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Abstract. We study a zero-range process with two species of interacting particles. 
We show that the steady state assumes a simple factorised form, provided the dynamics 
satisfy certain conditions, which we derive. The steady state exhibits a new mechanism 
of condensation transition wherein one species induces the condensation of the other. 
We study this mechanism for a specific choice of dynamics. 
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1. Introduction 

The zero-range process has recently received widespread attention in the literature of 
nonequilibrium stochastic processes [H [21 EH IH El El E] • There are three main reasons for 
this (see jS] for a review): (i) it provides a minimal model for more complex systems, 
(ii) the steady state has a simple factorised form, (iii) it exhibits a steady state phase 
transition, between a fluid and a condensate phase. The simple form of the steady state 
allows detailed analysis of this condensation transition. Apart from being of interest in 
its own right (as an example of a phase transition in an exactly solvable steady state), 
the transition is of particular interest because of its relation to coarsening phenomena 
in Id driven systems [21 H|- Further, it has provided a general criterion to determine the 
existence of phase separation in driven systems with conserved density [3] . 

In the zero-range process, particles hop from site to site with hop rates determined 
by the number of particles at the departure site. Two types of condensation transition, 
whereby a finite fraction of the particles condense onto a single site, have been studied 
jHj: condensation induced by site-dependent hop rates and condensation induced by 
the asymptotic dependence of the hop rate on particle number at the site. The former 
is related to Bose-Einstein condensation whereas the latter involves a spontaneously 
broken symmetry. 

Here, we generalise to a zero-range process with two species of particles. The 
hop rates of each species are determined by the numbers present of both species at 
the departure site. Under certain conditions which we specify, the steady state again 
assumes a simple factorised form, which we use to demonstrate a new mechanism of 
condensation transition. In this case the condensation of one species is induced by the 
distribution of particles of the other species. 

One perspective on this two species model is as a system of particles moving on an 
evolving landscape — the dynamics of one particle species depend locally on how many 
particles of the other species are present (the landscape) and these particles too are 
subject to locally determined dynamics. The new condensation mechanism then is one 
where the condensate is induced by the landscape whose evolution is in turn coupled 
to the particle distribution. Such interplay between particle dynamics and an evolving 
landscape is of considerable interest [HI E] • 

2. Steady State 

We define a two species zero-range process as follows. Consider a chain of L sites, with 
periodic boundary conditions, and containing N particles of species A and M particles 
of species B. Sites, labelled I = 1, . . . , L, contain ni particles of species A and mi 
particles of species B. Particles hop to their nearest neighbour site to the right, species 
A with rate u(n, m) and species B with rate v(n, m), where n and m are the occupation 
variables of the departure site. Thus the hopping rates, which depend on the number of 
particles of both species at the departure site, contain the interaction between species. 
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The steady state for the two species model can be derived in a similar way to that 
for the single species zero-range process [3]. First, we define P({ni}; {mi}) to be the 
steady state probability of finding the system in the configuration ({ni}; {mi}), where 
{n{} = rii,...,ni and {mi} = m ly . . . ,m L . Now, the steady state Master equation, 
which demands that the probability current into a particular configuration is equal to 
the probability current out of the same configuration, can be written 

L 

= ^2 + • • ,7^-1 + 1, raj-l, .,m { _i,mj, . . .) 

i=i 

-u{n h mi)P{. . . , nj_i, n h , mj_ ls m h . . .)} 9(ni) 

+ {u(n 1 _i,mj_ 1 + l)P(. . .,7ii_i,nj, ...;.. . , m f _i + l, m { -l, . . .) 

-v(nj, m/)P(. . . , nj_i, n h , m z , . . .)} 0(mj)] (1) 

where 0(x) is the usual Heaviside function and ensures that the site I is occupied in 
order for a particle either to have arrived or to be able to vacate there. The first term 
on the r.h.s of Q is a gain term due to an A particle hopping into site I from site I— 1, 
the second term is a loss term due to an A particle hopping out of site /; the third and 
fourth terms represent analogous processes for the B particles. We seek a solution to 
(HJ where P({n z }; {mi}) is given by a product measure (i.e. a product of factors) 

L 

P{{m}- {mi}) = z-^ M H f(m, mi) , (2) 

1=1 

where Z LjN ^ M is a normalisation. Equation ((TJ) can be satisfied by substituting in the 
steady state (J2J), and asking that the gain and loss terms due to the dynamics of the A 
particles cancel independently of the gain and loss terms due to the dynamics of the B 
particles. Then, if this cancellation is achieved for each term I in the sum individually, 
one finds that for all I, 

u(ni,mi)f(ni-i,mi-. 1 )f(ni,mi) = u(n^i+l, mj_i)/(n i _i+l, mj_i)/(nj-l, mi) , (3) 
for all ni 0, and 

v(ni,mi)f(ni-i,mi-i)f(ni,mi) = v(ni-i, m z _i+l)/(nz_i, mi-. 1 +l)f(n h mj-1) , (4) 

for all mi ^ 0. These equations in turn imply that 

u{n h mi)f(ni, m t ) u(n/_i+l, m/_i)/(nj_i+l, m/_i) 



f(m-l,mi) f(ni-i,mi-i) 
v(m,mi)f(m,mi) u(m_i,mi_i+l)/(nj_i,mi_i+l) 



constant , (5) 
constant , (6) 



f{n h mi-l) f(ni-i,mi-i) 
for all /. Both constants are set equal to unity (they only appear as an overall factor in 
the normalisation). The two relations © and (JBJ) imply a constraint on the choice of 
u(ni, mi) and v(jii, mi) i.e. we can use (J5J) and © to obtain two expressions for f(ni, mi) 
in terms of /(nj — 1, m/ — 1): 

,/ v /(n{-l,mj-l) /(ni-l,mi-l) , , 

/(n Is m z = — — = — — (7) 

u(ni,mi)v(ni-l,mi) u{n h mi-l)v{n h mi) 
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but both of these expressions must give the same result, therefore the hopping rates are 
required to obey the constraint, 

u(n h mj) v(n h mi) 
u{n u mi-l) v(ni — l,mi) 
for ni,mi 7^ 0. The choices of u(ni,0) and v(0,mi) remain unconstrained. 

Finally, the expression for f(ni,mi) is obtained by iterating (jSJ) and (JUJ), hence 

ni mi 

f(n hmi ) = ]\[u(i,m l )r 1 l[[v(0,j)r 1 , (9) 

i=\ 3=1 

where we have set /(0, 0) = 1. Although not immediately obvious, u(n,m) and v(n,m) 
do play symmetric roles in @, but this symmetry is obscured within the constraint (JBJ) 
on the hopping rates. 

Also note that we have the freedom to choose any desired form for f(ni,mi), then 
by substituting this form into (jSJ) and © we deduce the hopping rates which would lead 
to such a steady state — such rates are then guaranteed to satisfy the condition (jHJ). 

In order to look for phase transitions, it will prove useful to consider the 
normalisation, Z LtNjM , defined in equation (J2J). This quantity plays a role analogous to 
the canonical partition function of equilibrium statistical mechanics. It can be written 

L L L 

Z l ,n,m= S^ni-WS^mt-AfjJlffamt), (10) 

{ni},{mi} 1=1 1=1 1=1 

where the delta-functions ensure that we have N particles of species A and M particles 
of species B in the system. Using an integral representation of the delta-functions yields 

z L , N , M - 2. f^if^Ti z N + i yM+ i > (ii) 

{ni},{m t } 

which leads us naturally to define the generating function 

N M 

F(^,y) = ^^^ m /(n,m). (12) 

n=0 m=0 

We now assume that in the limit N, M, L — > 00 (with the particle densities held 
fixed), the integral in equation (fTTjl is dominated by its saddle point. The saddle point 
equations, which involve the particle densities of species A and B, pa = N/L and 
Pb = M/L, are found to be 

d d 
Pa = z—\nF(z, y) , p B = y—lnF(z, y) . (13) 

These equations determine z and y. For the saddle point to be valid, z and y must be 
less than or equal to the radii of convergence of F(z, y) in order that we can perform 
the sum (112)1 in the first place. We note also that because all derivatives of F(z,y) are 
positive, the saddle point, if valid, must be unique. In the following, we will find that 
it is only possible to solve (fTB*)) in the allowed ranges of z and y for certain values of pA 
and pb- 
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3. Phase Transition — analysis 

We now consider how a condensation transition may arise from the interaction of the 
two species. A general case of interest is where the hopping rates of one of the particle 
species depend only on the number of the particles of the other species at the site e.g. 
v(n,m) = r(n). Then, by Q, f(n,m) = r(n)~ m s(n), and u(n,m) is determined by (JSJ) 
(r(n) and s(n) are general functions of n). In this case, the A particles play the role of 
an evolving landscape that determines the dynamics of the B particles. For simplicity, 
we consider the case where r(n) — 1 + c/(n + 1) and s(n) = 1, then 

v(n, m) = 1 + ^q-y , u(n,m) = ^ - + " £ J , (14) 

for all values of n and m; c is a constant. Physically, these hopping rates are such 
that the B particles hop more slowly as the number of A particles at the departure 
site increases, while the A particles slow as the number of B particles increases but 
they hop more quickly the more particles of the same species occupy the departure site. 
The generating function and its first derivatives, having performed the sum over m, are 
obtained as 

oo 1 _|_ _|_ 

' {l+n)(l — y) + c 



n=0 



z^-F(z, y) = Y nz n - l t C + U , , ( 16 ) 

dz 1 ,y) ^ (i + „)(i_ y ) + c ' ^ > 



dy ^ [(1 + 0(1 -y) + c} 2 

The radii of convergence of these series are z = 1 and y = 1 . These expressions determine 
z and y, given and as prescribed by the saddle point equations 

One can demonstrate the existence of a phase transition in the following way. For a 
given z, both p^ and are monotonically increasing functions of y. Then by analysing 
the saddle point equations for the four cases y = 0, z = 0, y = 1, and z = 1, the 
dependences of pa and pB on y, for fixed z, are obtained. These are shown in figure 
where < z\ < < 1. Now, if we consider a system containing a density of A 
particles, then the solution of (fT^j) requires that z and y lie in the range z\ < z < Z2 and 
< y < 1. In this range, increases monotonically from p# = 0, where z = Z2 and 
y = 0, to its maximum value allowed by the saddle point equations, where y = 1 and 
z = 2Tx. Therefore we can only solve the saddle point equations for p^ and ps if Pb is 
less than this (p^-dependent) maximum. When p# exceeds this maximum, the saddle 
point approximation is no longer valid and a phase transition ensues. The critical line 
is given by the values of pa and p#, as a function of z, where y — 1. Along this line pa 
and p# are related by 

p B = (l + p A )/c. (18) 
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Figure 1. Schematic dependences of the particle densities pa and ps for contours of 
fixed z and as a function of y. The dashed line in the right hand graph illustrates how 
Pb varies as a function of z and y given that pa is fixed. 



Pa 




1/c Pb 

Figure 2. Phase diagram for hopping rates 114|) . 

This is most easily seen by setting y — 1 in (fTB]). (JTHJ) and (|T7j). and then spotting that 
one can write ydF(z,y)/dy = [F(z,y) + zdF(z,y)/dz]/c. Then dividing through by 
F(z, y) yields (|18|) . Hence we obtain the phase diagram shown in figure ED 

We interpret this as a transition between a fluid phase, where the saddle point 
equations can be solved, and a condensate phase. In the condensate phase, ps exceeds 
a maximum critical value and the excess B particles condense onto a single site. The 
number of B particles in the condensate is proportional to the system size L — in the 
limit N, M, L — > oo it contains a finite fraction of an infinite number of particles. 
However, we note that for m — > oo, u(n,m) — > if n is finite. But the current of 
A particles must be the same across all bonds in the system as a result of continuity, 
and this current must be finite — the current of A particles away from the condensate is 
finite — therefore we must have that n —>■ oo at the condensate in order that the current 
out of the condensate site is finite. For n large, u(n,m) ~ exp(m/n 2 ) and this must be 
finite, therefore m ~ n 2 at the condensate then since m oc L here, n oc L 1 / 2 . Thus the 
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P(n) P(m) 




Figure 3. Log-log plot of on the left P(n) vs. n, and on the right P(m) vs. m, for 
systems of size L — 100 and c = 2. The circles correspond to densities pa — 1/2 = pb', 
the crosses correspond to densities pa = 1/2 and ps = 7. 

condensate phase is characterised by a condensate of £> particles which is sustained by 
a 'weak' condensate of A particles. This picture is consistent with numerical results: 
an exact recursion relation for the partition function can be obtained from ()10|) and is 
found to be 

N M 

Z L ,N,M = ^2 /L, m ) Z L-l,N-n,M-m ■ (19) 

n=0 m=0 

This expression was iterated exactly on a computer up to systems of size L = 100, for 
different values of the densities and using c = 2. It is also straightforward to obtain an 
exact expression for P(n,m), the probability that a site contains exactly n particles of 
species A and m particles of species B. P(n, m) is given by 

P(n, m) = Z^ NM Z L _ ltN _ ntM _ m f(n, m) . (20) 

By summing this equation over norm one obtains equations for P(n), the probability 
of finding exactly n particles of species A at a site, and P(m), the probability of 
finding exactly m particles of species B at a site. These probabilities are plotted in 
figure 01 for two different values of the densities. The circles correspond to densities 
Pa = 1/2 = pb and the system is in the fluid phase — both P{n) and P(m), as expected, 
decay exponentially. For pa = 1/2 and ps = 7, represented by the crosses, the system is 
in the condensate phase. Within finite size effects, these results are consistent with the 
following picture: the condensate of B particles (containing ~ L particles) is sustained 
by a 'weak' condensate of A particles (containing ~ L 1 ^ 2 particles); the condensates 
exist on a power-law distributed background of B particles forming a critical fluid, and 
an exponentially distributed background of A particles forming a normal fluid. 
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4. Conclusion 

We have shown that one may define a two species zero-range process where the steady 
state is described by a product measure (J2J) and we have derived conditions under which 
the product measure holds (jHJ). This simple form for the steady state makes the model 
amenable to detailed analysis. In particular, we have shown that the model undergoes 
a novel kind of condensation transition. 

It is clearly of interest to investigate the generality of this transition with respect to 
the hopping rates. Since we can always construct dynamics, using (0121), that lead to a 
product measure (J2J) with arbitrary f(n,m) there is scope within the model for further 
mechanisms of condensation transition ^U]. Further analysis of the nature of the weak 
condensate is also desirable. 
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